133 



Landscape Ecology 17 : 133-155, 2002. 

© 2002 Kluxver Academic Publishers. Printed in the Netherlands. 


Recent forest cover type transitions and landscape structural changes in 
northeast Minnesota, USA 

Peter T. Wolter 1 ’ 2 ’* and Mark A. White 1 

1 Center for Water and the Environment, Natural Resources Research Institute, University of Minnesota 
Duluth, Duluth, MN 55811, USA; 2 Current address: University of Wisconsin - Green Bay, 2420 Nicolet Drive, 
Green Bay, WI 54311, USA; * Author for correspondence (e-mail: Pwolter@nrri.umn.edu) 

Received 1 January 2001; accepted in revised form 7 September 2001 


Key words: Change detection, Forest classification, Forest dynamics, Forest fragmentation, Landsat, Landscape 
indices, Landscape pattern 


Abstract 

Landsat TM satellite data covering an approximate 5-year interval (1990-1995) were used to quantify spatial 
pattern and transition rates between forest ecological states for a 2.76 million ha region in northeast Minnesota. 
Changes in forest cover were stratified by Ecological Subsection, management status, and by ownership catego¬ 
ries using a 1995 digital ownership layer. Approximately 4.2% of the 1990 mature forested area was converted to 
early successional types by 1995. Of this 4.2%, private lands accounted for 33%, federal lands 31%, county 
lands 20% and state lands 16%. Notable conversion percentages by cover type category were spruce-fir (-5.3%), 
aspen-birch (-4.7%), jack pine (-4.6%) and black spruce (-3.0%). Transition rates were also adjusted to fit ten- 
year time intervals. Shannon-Weaver Eveness and edge density of cover types increased over the study period as 
relative contagion and interior forest area decreased. These trends suggest both smaller patches and a more even 
distribution of cover types. Area of upland conifers, lowland conifers and lowland hardwoods decreased while 
the area of mature upland hardwoods increased in most patch size classes except the > 500 ha class which showed 
a substantial decrease in area. The area of early successional types increased in most patch size classes. Non¬ 
industrial private forestland had the lowest proportion of interior forest of all ownership categories - decreasing 
by 13.5% in five years. Smaller average cut-unit sizes and uncoordinated forest management is the likely cause 
since cutting rates between private and public forestland were similar. 


Introduction 

Minnesota is expected to increase timber harvesting 
by roughly 25% in coming decades to fuel the fiber 
demands of the paper industry (Jaakko 1994). Unfor¬ 
tunately, one of the problems in monitoring landscape 
change in this region is that the current rate of forest 
cutting is faster than the forests can be re-inventoried 
on the ground (Czaplewski 1999). For example, the 
ground-based U.S. Forest Service’s Forest Inventory 
and Analysis (FIA) program is repeated only once 
every 10-15 years (Czaplewski 1999). Consequently, 
over one half of the latest FIA database is out of date 
at any one time, and FIA data greater than five years 
old are, therefore, not regarded as reliable (American 


Forest Council 1992). Even if the FIA database was 
current, these data are made up of point samples and 
are not well suited for studying spatially explicit for¬ 
est change dynamics on a landscape basis. 

One solution to the inventory problem is the use 
of space-borne sensors such as Landsat Thematic 
Mapper (TM). TM data have been found to be an ef¬ 
ficient tool for the detection of forest canopy changes 
due to logging operations (Olsson 1994), insect infes¬ 
tation (Vogelmann and Rock 1989), acid rain (Vogel- 
mann and Rock 1988) and forest succession (Hall et 
al. 1991). Our objectives were to determine if we 
could detect changes in forested landscape structure 
using Landsat TM data. We hypothesized that 
changes in forest spatial structure over time will vary 
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Figure 1. Study area in northeast Minnesota, two Landsat TM footprints that cover the area, and Ecological Subsection boundaries. 


by ownership, forest management status, and ecolog¬ 
ical subsection. 

Study area 

The 2.76 million hectare study region in northeastern 
Minnesota is composed largely of boreal forests and 
includes Voyageurs National Park, Superior National 
Forest, and the protected wilderness of the Boundary 
Waters Canoe Area (BWCA) (Figure 1). Forest man¬ 
agement practices within this region range from in¬ 
tensively managed private industrial forests to the un¬ 
managed BWCA and Voyageurs National Park. It 
should be noted that logging has not been permitted 
in the BWCA since 1978. The BWCA historical dis¬ 
turbance regime, driven mainly by large fires (Hein- 
selman 1973, 1996) and wind throw (Frelich and Car¬ 
lton 2000), is vastly different from the disturbance 
regime caused by timber harvesting outside wilder¬ 
ness areas (Hall et al. 1991). Heinselman (1973); Hall 
et al. (1991); Pastor and Mladenoff (1992); Pastor et 
al. (1993); Frelich and Reich (1995) have previously 
described the ecology of this study region. 


Satellite data 

Two satellite footprints were needed to cover the 
study area: Worldwide Reference System (WRS) 
path/row 27/27 and 26/27 (Figure 1). TM data from 
12 May 1990, 31 July 1990, 24 February 1991, 20 
September 1991 and 29 July 1995 were used for path/ 
row 27/27. Path/row 26/27 TM data included 4 March 
1984, 26 May 1986, 2 July 1988, 26 September 1990 
and 1 June 1994. 

Summer TM data from 31 July 1990 (WRS 27/27) 
and 2 July 1988 (WRS 26/27) (hereafter nominally 
referred to as 1990) were calibrated to reflectance ac¬ 
cording to Price (1987) and geometrically registered 
to UTM zone 15 (NAD83) coordinates. We achieved 
a root-mean-squared (RMS) error of unit weight of 
approximately 0.37 pixels (28.5 meter pixels) for the 
fit between digital Landsat TM data and the 1:24,000 
scale USGS topographic maps using 39 evenly dis¬ 
tributed ground control points. 

Each of the remaining scenes was calibrated to re¬ 
flectance then radiometrically rectified (Hall et al. 
1991) to the July TM data for each respective foot¬ 
print. These data were then coregistered to the July 
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data. For better spatial continuity of the TM data 
through time, a coregistration algorithm which as¬ 
signs output pixel locations based on interpolating in¬ 
put brightness values in two orthogonal directions (bi¬ 
linear interpolation) was used. Nearest neighbor im¬ 
age registration algorithms preserve precise input 
pixel spectral information but tend to corrupt linear 
features due to poor spatial relocation of pixels, 
which may contribute a significant amount of error in 
multi-temporal studies (Logan and Strahler 1979). 
Assessments of the image-to-image registrations were 
performed independently by selecting 26 checkpoints 
per scene. All RMS errors were within 0.5 pixels (first 
order) for the fit between the July dates from 1990 
and all other Landsat data sets. 

Forest classification 

The classification techniques employed in this study 
are essentially those used by Wolter et al. (1995) in 
northern Wisconsin, USA. Therefore, only a brief de¬ 
scription of forest classification methods will be dis¬ 
cussed. Forest canopies in this region are phenologi- 
cally similar during mid-summer, and hence, satellite- 
based classification aids gleaned from winter, spring, 
and autumn phenology were made via comparison to 
July TM data. 

July TM images from each footprint were classi¬ 
fied into conifer, hardwood, mixed conifer-hardwood, 
and non-forest using Minnesota Department of Natu¬ 
ral Resources (MnDNR) Phase II forest inventory 
data (Minnesota Department of Natural Resources 
1995) and 1:40,000 color infrared (CIR) aerial pho¬ 
tographs (May and September 1991) as ground truth. 
The hardwoods class was used to clip the September 
and May Landsat data for each footprint to classify 
black ash ( Fraxinus nigra), maple-dominated ( Acer 
saccharum and A. rubra) northern hardwoods and as¬ 
pen-birch ( Populus tremuloides-Betula papyrifera ) 
cover types. The conifers class was used to clip the 
winter TM data for classifying tamarack ( Larix 
laracina ), a deciduous conifer. Winter TM data were 
also used to detect the presence of understory conifer 
components among types classified as hardwood us¬ 
ing summer TM data. 

Peak maple senescence in this region typically oc¬ 
curs near the autumnal equinox (Eder 1989) while 
black ash trees are completely defoliated by this time 
(Wolter et al. 1995). Conversely, aspen-birch is the 
first hardwood cover type to flush leaves in the spring 


- usually by mid-May (Ahlgren 1957). Senescent, 
maple-dominated, northern hardwood cover types 
were classified using a combination of July and Sep¬ 
tember TM data. Black ash, aspen-birch and tamarack 
types were classified by applying thresholds to nor¬ 
malized difference vegetation index (NDVI) differ¬ 
ence images: September-July, May-July and winter- 
July respectively. NDVI is defined as: 

((TM4 - TM3)/(TM4 + TM3) + 1) x 100 (1) 

MnDNR Phase II forest inventory and CIR aerial 
photos were used to guide threshold selection for the 
difference images. Remaining cover types were clas¬ 
sified using traditional iterative classification tech¬ 
niques using Phase II forest inventory, CIR aerial 
photos and field visits for validation (see Wolter et al. 
(1995)). 

Overall classification accuracy for the WRS 27/27 
footprint was 75% for 37 classes (KHAT 0.74) (Ap¬ 
pendix 1). The error matrix for the WRS 26/27 foot¬ 
print (not presented) was similar to the 27/27 foot¬ 
print with an overall accuracy of 73% (KHAT 0.72). 
Accuracy assessment methodologies for both foot¬ 
prints were performed according to Wolter et al. 
(1995) where MnDNR Phase II forest inventory and 
1995 CIR aerial photos were used as ground truth. 

Change detection 

TM scenes 7/29/1995 and 6/1/1994 (hereafter nomi¬ 
nally referred to as 1995) represented the second 
time-step for footprints WRS 27/27 and 26/27 respec¬ 
tively (Figure 1). To detect land-cover changes in 
1995, we first identified pixels that exhibited a differ¬ 
ence, either positive or negative, in reflectance be¬ 
tween 1990 and 1995. To accomplish this, we utilized 
a normalized version of the short-wave infrared/near 
infrared ratio (SWIR/NIR) using TM bands centered 
at 1.65 jiim (TM Band 5; 1.55-1.75 pm) and 0.83 jitm 
(TM Band 4; 0.76-0.90 pm) respectively. The nor¬ 
malized difference SWIR/NIR ratio (NDSN) is de¬ 
fined as: 

((TM5 - TM4)/(TM5 + TM4) + 1) x 100 (2) 

This basic ratio has been used extensively in the east¬ 
ern United States to detect changes in forest cover due 
to insect outbreaks (Vogelmann and Rock 1989) and 
forest decline (Vogelmann and Rock 1988). Vogel- 
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mann (1990) showed that the SWIR/NIR ratio is far 
more sensitive to forest canopy disturbance than the 
widely used NDVI. Interpreted CIR aerial photo¬ 
graphs (May 1991 and September 1997) were use to 
determine appropriate discriminant values for change 
in the NDSN difference images. 

After change pixels between 1990 and 1995 were 
identified, they were used as a template to clip the 
classification from 1990. The clipped 1990 classifica¬ 
tion was then used as a template to stratify and clas¬ 
sify raw TM data from 1995 on a class-by-class ba¬ 
sis. For example, change pixels originating from the 
1990 aspen-birch class were used to focus the classi¬ 
fication of raw TM data for 1995. In doing this, po¬ 
tential illogical transition classes for this time span 
may be ruled out in favor of more logical aspen-birch 
transition endpoints (e.g., upland grass/brush, regen¬ 
eration, burned or flooded) when 1995 spectra are 
confused with illogical transition types. An aspen- 
birch stand that was harvested between 1990 and 
1995 for example, should not be classified as erica- 
cious brush, sphagnum moss or acid-bog conifers in 
1995. This would be a highly unlikely transition of 
cover types within the five-year time frame. Con¬ 
versely, if some of the change identified by the dif¬ 
ference image classifies to mature hardwood it is most 
likely aspen-birch and can be treated as no change. 
Once all cover type transitions were identified, they 
were merged with the 1990 classification to produce 
a complete classification layer for 1995. 

Accuracy assessment for 1995 change classes (Ap¬ 
pendix 2) was accomplished by randomly selecting 
349 reference sites from interpreted CIR aerial pho¬ 
tographs from September 1997. Overall accuracy of 
the change classes was 89% (KHAT 0.87). 

The spectral properties of harvested forests, and 
forest canopies damaged by insects can be difficult to 
discern using traditional classification algorithms. 
Qualitatively though, distinguishing between har¬ 
vested forests and natural forest disturbance is rela¬ 
tively simple. Harvesting creates distinct geometric 
shapes with sharp gradients between forest and non¬ 
forest in northern Minnesota where clearcutting is the 
primary mode of wood removal (Plate 1). Soft, con¬ 
tinuous gradients or highly irregular boundaries of 
defoliation/deforestation characterize natural disturb¬ 
ance. However, much of the timber that was blown 
down between 1990 and 1995 in the northeastern por¬ 
tion of this study area was harvested by salvage op¬ 
erations (Wayne Russ, USFS, personal comm.). 


Landscape analysis 

Landscape stratification 

Since environment (soil, landform and climate), own¬ 
ership and timber management status can contribute 
to spatial heterogeneity (see Crow et al. (1999)), we 
analyzed landscape changes by stratifying our analy¬ 
sis by ecological land units in the MN DNR Ecologi¬ 
cal Classification (Minnesota Department of Natural 
Resources 1999) and by ownership (Gap Stewardship 
Database, Minnesota Department of Natural Re¬ 
sources (1998)). Ownership categories were stratified 
by timber management status (managed, unmanaged 
and unknown). In our study, the term “unmanaged” 
refers to timber or forested areas in reserves, such as 
the BWCA, parks and natural areas, which restrict 
traditional forest harvesting operations. “Unknown” 
timber management status lands are non-industrial 
private forests. 

We analyzed the study area overall and by subsec¬ 
tion (Figure 1). Detailed analyses were focused on 
one subsection, the Border Lakes, by ownership and 
timber management status. This is the only subsection 
that contains significant amounts of mature forest 
cover under each timber management status category. 
Other subsections contain less than 5% of land that 
can be identified as unmanaged (Gap Stewardship 
Database, Minnesota Department of Natural Re¬ 
sources (1998)). Furthermore, the large area of wil¬ 
derness (BWCA) with a natural disturbance regime 
within the Border Lakes Subsection serves as a land¬ 
scape control for this subsection. 

Transition matrix analysis 

To quantify vegetation transitions we utilized the GIS 
overlay procedures in ERDAS Imagine 8.4 (ERDAS 
1997). Cover types from a given state (time t) are 
used to clip out cover types from the next state (time 
t + 1). From these data we generated matrices of 
transition probabilities. Methods for calculating tran¬ 
sition probabilities follow Pastor et al. (1993). The 
transition probabilities for time interval t were calcu¬ 
lated by: 

m 

Pij, t = n i jX, n ij ( 3 ) 

J = 1 

where p ijt is the probability that a given area of the 
landscape has changed from class i to class j during 
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Plate 1. A. 1988 TM false color composite of a portion of the Superior National Forest and the Pat Bayle State Forest in northeast Minne¬ 
sota. Image was produced using TM 1.65, TM 0.83 and TM 0.56 /xm bands in the order of red, green and blue. 1 = Grand Marais, MN, 2 
= Greenwood Lake and 3 = Boundary Water Canoe Area Wilderness. B. False color composite of the same area from 1994 - pink areas are 
forest cuts. C. 1988 to 1994 change image produced using 1994 TM 1.65, 1988 TM 1.65 and 1988 TM 0.56 /xm bands in the order of red, 
green and blue. White areas are forest cuts from 1988 and earlier. Red areas are forest cuts that occurred between 1988 and 1994. 


time interval t. n tj is the land area undergoing the 
transition from i to j across the landscape of m classes 
or states. In this analysis, the sum of p tj = 1, indicat¬ 
ing class / moves with some probability value to an¬ 
other class, including itself (Johnson and Fryer 1987). 
Transition probabilities over the time interval t were 
normalized to a 10-year interval for the periods from 
1988-1994 and 1990-1995 using the equation de¬ 
scribed by Pastor et al. (1993): 

Pj j = 1 — e (ln(1 ~ 7r))// when i A j (4) 

m 

P>,i= ] ~ 'Zp.j ( 5 ) 

J = 1 

Equation (4) applies to off-diagonal elements of the 
transition matrix. Equation (5) applies to cases when 
i = j , or on-diagonal elements and represents reten¬ 
tion probabilities. 

Landscape metrics 

The software program APACK was used to calculate 
the following landscape statistics (Mladenoff and De- 


Zonia 1999). To examine overall changes in land¬ 
scape structure and pattern, we used two indices 
based on information theory. Shannon-Weaver Even¬ 
ness (SWE) measures the evenness of the proportional 
distribution of patch type area. SWE increases with 
increasing evenness (McGarigal and Marks 1995). 

H' 

SWE = - (6) 

maximum diversity possible 

Where measured diversity (//') = 

classes present 

- X P(i) x ln(p(/)) (7) 

i = 1 

and p equals the proportion of patch type i and maxi¬ 
mum possible diversity equals In (<classes present). 
This measure describes landscape structure, but has 
no spatial component (Li and Reynolds 1994). 

Relative contagion ( RC) was used to examine the 
spatial pattern or texture of the study landscape. This 
index is derived from the cover type co-occurrence 
probabilities of the cover type adjacency matrix, 
which quantifies the frequency that two cover types 
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share a common border. RC measures the evenness 
of the adjacency matrix; it quantifies the degree of 
clumping or aggregation of cover types on the land¬ 
scape. Values range between 0 and 1; high values in¬ 
dicate with a high degree of aggregation with fewer, 
large patches, while low values indicate landscapes 
composed of smaller patches (Li and Reynolds 1993). 

RC=1.0-( -—- \ (8) 

\maximum possible diversity ) 

n n 

W'=-H pMp>) (9) 

i= lj= 1 

Where n is the total number of classes in a landscape, 
and p x j is the probability of patch type i being adja¬ 
cent to patch type j. Maximum possible diversity 
equals: 

2 x ln(classes present) (10) 

To analyze finer scale changes in forest patch 
structure, we calculated a set of patch statistics based 
on patch size and size class distribution. Patch size 
class distribution analysis has proven useful in ana¬ 
lyzing changes and differences in forested landscape 
pattern and structure (Mladenoff et al. 1993). Vegeta¬ 
tion patch types were aggregated into six classes: up¬ 
land conifer (UC), lowland conifer (LC), upland hard¬ 
wood (UH), lowland hardwood (LH), regeneration 
(REGEN) and upland grass-brush (UG-UB). The fol¬ 
lowing statistics were calculated for each class for 
each time period for the landscape stratification de¬ 
scribed above: number of patches, mean patch size 
(ha) and standard deviation, and largest patch. We 
analyzed patch size class distribution by calculating 
the number of patches and the area for the following 
size classes in hectares: 0-1, 1-5, 5-10, 10-50, 50- 
100, 100-500 and > 500. 

Mean edge-to-edge patch nearest neighbor (MNN) 
statistics were calculated on early successional forest 
patches to analyze spacing differences by ownership 
category. Ten stratified random subsamples per cate¬ 
gory (mean subsample size ~ 10,000 ha) were se¬ 
lected to avoid fragmented ownership patterns. A one¬ 
way analysis of variance (ANOVA) on log trans¬ 
formed MNN data was performed to determine if 
mean patch spacing distances from at least one of the 
ownership categories was different from the rest. If 
so, pairwise multiple comparisons using the Tukey 


Test would be used to assess patch spacing differ¬ 
ences between ownership categories. 

Edge Density (ED) measures the edge length per 
unit area in a given landscape. The ED function that 
APACK uses corrects for pixel shape and saw-tooth 
polygon edges that are characteristic of raster data. 
We measured edge density for two dates for the six 
cover classes and total edge density for seven subsec¬ 
tions and three management classes (managed, un¬ 
known, unmanaged) for the Border Lakes Subsection. 
We calculated percent difference from 1990 to 1995. 

Interior Forest Area (IFA) is used as an indication 
of the extent of large forest patches in a landscape 
(Sachs et al. 1998). IFA is measured by determining 
the amount of forest area that is beyond a certain dis¬ 
tance from a non-forest edge (e.g., grass, brush, roads, 
water). We aggregated all forest types, except regen¬ 
erating forest classes, into one class then calculated 
IFA, according to Sachs et al. (1998), using a 100-m 
forest edge buffer. We calculated IFA overall and by 
ownership and management status categories for 
1990 and 1995. We also calculated the proportion of 
IFA to total forest area per time step. 

Results 

Cover type classes from 1990 and 1995 were recoded 
from 37 classes to 23 classes (see Appendices 1 and 
3) prior to calculation of transition matrices. The ag¬ 
gregation of classes both simplified analyses and el¬ 
evated overall accuracy figures to 79% (KHAT 0.772) 
and 77% (KHAT 0.749) for WRS paths 27/27 and 
26/27 respectively. Then 5-year (not presented) and 
10-year transition matrices (Appendix 3) were created 
using output from the MATRIX program (ERDAS 
1997). 

Five-year transitions 

Of all 1990 mature forest, trees pole size and larger, 
cover types combined (1.76 million ha), 1.9% of the 
area became native grasses, 1.8% transitioned into 
hardwood regeneration and 0.1% went into conifer 
regeneration by 1995. Transitions to flooded, assumed 
to be naturally induced changes, were minimal. Of the 
forested area in 1990, 4.2% was harvested, flooded or 
burned by approximately 5 years later (Plate 1, Ta¬ 
ble 1). The percent-harvested value was calculated by 
subtracting cover type transitions to flooded, emer¬ 
gent, water and burned from the total (Table 2). Com- 
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Table 1. Mature forest area changed to early successional types by ownership and timber management status category between 1990 and 
1995 for the whole study area. 




1990 mature forest 

area (ha) 

Forest changed by 

1995 (ha) 

1990-1995 forest 

change (%) 

Annual rate (%) 

Federal 

managed 

422,894 

20,180 

4.77 

0.95 


unmanaged 

307,817 

2,647 

0.86 

0.17 

State 

managed 

195,202 

10,598 

5.43 

1.09 


unmanaged 

57,549 

1,029 

1.79 

0.36 

County 

managed 

289,489 

14,835 

5.12 

1.02 


unmanaged 

3,148 

34 

1.08 

0.22 

Private Industrial 

managed 

37,764 

3,139 

8.31 

1.66 


unmanaged 

9,948 

1,006 

10.11 

2.02 

Non-indust. Private 

unknown 

415,059 

19,721 

4.75 

0.95 

Tribal 

managed 

16,499 

456 

2.76 

0.55 

All categories combined 

1,755,367 

73,645 

4.20 

0.84 


Table 2. Percent change in mature forest area for the whole region 
from 1990 to 1995 by cover type and percentage of that change 
due to forest harvesting operations. 



% change 

% of change harvested 

jack pine 

-4.55 

92.0 

red pine 

-2.93 

96.0 

spruce-fir 

-5.33 

95.2 

cedar 

-2.34 

88.2 

tamarack 

-2.31 

81.3 

black spruce 

-3.04 

88.1 

bog conifers 

-1.96 

83.1 

black ash 

-3.19 

75.7 

aspen-birch 

-4.72 

97.6 

northern hwd. 

-2.72 

99.4 


bining human and natural disturbance, notable 
changes included a 5.3% decrease in the spruce-fir 
cover type, 4.7% decrease in the aspen-birch cover 
type, 4.6% decrease in the jack pine cover type and a 
3.0% decrease in the black spruce cover type. 

In 1990, federal lands made up the greatest pro¬ 
portion of mature forest area in the study area fol¬ 
lowed by private, county and state lands (Figure 2a). 
Of the managed federal forestlands, transitions from 
mature forest to non-forest and early successional for¬ 
est accounted for 27.4% of the total forest disturbed 
(i.e., cut, flooded or burned) in the study area by 1995 
(Figure 2b). State lands contributed 14.4% of the to¬ 
tal and county lands 20.1% of the total mature forest 
cover disturbed by 1995. By 1995, non-industrial pri¬ 
vate lands accounted for 26.8% of the total decrease 
in mature forest area (Figure 2). Managed state lands, 


with 11.1% of the total forest area, accounted for 
14.4% of the total loss of mature forest area. 

Table 1 shows a breakdown by ownership and 
management status categories of mature forest area 
converted to early successional types between 1990 
and 1995. Mature, federally managed forest de¬ 
creased 4.8% in area. Managed state and county ma¬ 
ture forest decreased 5.4% and 5.1% in area, respec¬ 
tively. By 1995, 4.8% of the mature, non-industrial 
private forest was disturbed due to logging, wind 
throw, flooding or burning. The area of mature pri¬ 
vate industrial forests, although a small percentage of 
the total forested area in this region, decreased by 
8.3%. The ~ 10,000 ha of unmanaged, private indus¬ 
trial forest is primarily made up of lands dedicated to 
mining operations. Nevertheless, mature forest on 
these lands decreased by 10.1% in area by 1995. 

Landscape pattern and structure 

Entire study area 

Shannon-Weaver Evenness (SWE) increased for all 
subsections from 1990 to 1995. Increases ranged from 
2.3 to 4.5% (Table 3). Increasing SWE values indi¬ 
cate that the distribution of area by cover type is be¬ 
coming more even through time. Relative contagion 
showed a decrease, with values ranging from -5.4 to 
-11% (Table 3). An overall decrease in RC suggests 
both smaller patches and higher evenness of cover 
type distribution. The North Shore Uplands subsec¬ 
tion (Figure 1) had both the smallest SWE and larg¬ 
est RC values of all the subsections studied, while 
Littlefork-Vermillion Uplands had the highest SWE 
and lowest RC values. The Border Lakes, Toimi Up- 
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Figure 2. A. Breakdown of 1990 mature forest area (1,755,367 ha) by ownership and management status. B. Percent of mature forest area 
that was converted to early successional types (73,645 ha) by 1995. 


lands and Laurentian Uplands subsections were all 
intermediate with 1995 SWE and RC values —0.81 
and 0.27 respectively. Based on SWE and RC, there 
appears to be three distinct groups of subsections: 1) 
North Shore Uplands, 2) Border Lakes-Laurentian 
Uplands-Toimi Uplands and 3) Nashwauk Uplands- 
Tamarack Lowlands-Littlefork-Vermillion Uplands 
(Table 3). 

Border lakes subsection 

Within the Border Lakes Subsection, management 
status showed a strong effect on SWE and RC similar 
to the overall study area. Managed and unknown 
showed similar changes with SWE increasing by 4.4 
and 3.4% respectively while unmanaged showed a 
relatively small increase of 0.60%. RC values reveal 
a similar trend with changes of -11.0 and -9.3% on 
managed and unknown, while unmanaged decreased 
by 1.4% (Table 3). The BWCA wilderness contains 
— 91 % of the unmanaged mature forest area in this 
ecoregion. 


greatly (from 24 to 58% for regenerating forest, and 
from 9 to 17% for upland grass-brush). Upland hard¬ 
wood forest showed the least amount of change from 
1990 to 1995 (Table 4). The North Shore Uplands 
subsection had the lowest overall edge density (1995) 
of all subsections, while Tamarack Lowlands and 
Nashwauk Uplands were the highest. Overall subsec¬ 
tion edge densities point to the same three distinct 
groups as SWE and RC (Table 4). 

Border lakes subsection 

Managed and unknown categories within the Border 
Lakes Subsection showed edge density changes simi¬ 
lar to the overall study area with upland and lowland 
conifer and lowland hardwood decreasing and early 
successional types increasing (Table 5). Edge densi¬ 
ties in unmanaged land, which is largely wilderness 
(BWCA), changed little across classes with the ex¬ 
ception of regenerating forest, which increased by 
42%, however, regenerating forest covers only 0.4% 
of the unmanaged landscape (Table 5). 


Edge density 


Cover type area by patch size class 


Entire study area 

Edge densities of upland and lowland conifer, and 
lowland hardwood forest decreased slightly but con¬ 
sistently while those of early successional types (re¬ 
generating forest and upland grass-brush) increased 


Entire study area 

The area of upland conifer, lowland conifer, and low¬ 
land hardwood forest decreased for all patch size 
classes, with the exception of a small increase in the 
area of 0-1 ha in upland conifer forest (Figure 3). 
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Table 3. A) Shannon-Weaver Evenness (SWE) and B) Relative 
Contagion (RC) and percent difference (% diff) for 7 subsections 
for the two time periods. C) SWE and D) RC and percent differ¬ 
ence for three timber management classes and two time periods for 
the Border Lakes Subsection. 


A) 

SWE 

1990 

1995 

% diff 

Border Lakes 

0.758 

0.786 

3.69 

Toimi Uplands 

0.785 

0.820 

4.46 

Laurentian Uplands 

0.786 

0.816 

3.82 

North Shore Uplands 

0.688 

0.710 

3.17 

Nashwauk Uplands 

0.840 

0.866 

3.10 

Tamarack Lowlands 

0.855 

0.875 

2.34 

Littlefork-Vermillion 

0.870 

0.897 

3.10 

B) 




RC 

1990 

1995 

% diff 

Border Lakes 

0.300 

0.277 

-7.67 

Toimi Uplands 

0.291 

0.262 

-9.97 

Laurentian Uplands 

0.294 

0.268 

-8.84 

North Shore Uplands 

0.355 

0.336 

-5.35 

Nashwauk Uplands 

0.251 

0.225 

-10.36 

Tamarack Lowlands 

0.255 

0.233 

-8.63 

Littlefork-Vermillion 

0.245 

0.218 

-11.02 

C) 




SWE 

1990 

1995 

% diff 

unmanaged 

0.664 

0.668 

0.60 

unknown 

0.800 

0.828 

3.50 

managed 

0.801 

0.836 

4.37 

D) 




RC 

1990 

1995 

% diff 

unmanaged 

0.361 

0.356 

-1.39 

unknown 

0.270 

0.245 

-9.26 

managed 

0.273 

0.243 

-10.99 


Upland conifer forest showed the greatest decrease 
(8%) in the largest size class (> 500 ha). Thirty four 
percent of the upland conifer area was in the >500 
ha class in 1990. Upland hardwood forested area in¬ 
creased in all classes except the largest, which de¬ 
creased by 13% (Figure 3a). The area of forest regen¬ 
eration increased in all size classes with values rang¬ 
ing from 20% in the 0-1 ha class to 125% in the 50- 
100 ha class. The area of upland grass-brush 
increased by ~ 10% in the first 3 size classes, de¬ 
creased in the 10-50 and 100-500 classes and in¬ 
creased by 45% in the > 500 ha class (Figure 3b). 


These data show a general trend of decreasing area 
in most patch size classes for upland conifer, lowland 
conifer and lowland hardwood, while upland hard¬ 
wood increased in most size classes except the > 500 
ha class which showed a substantial decrease in area. 
The early successional types showed corresponding 
increases in most patch size classes. 

Ten-year transition probabilities for the study area 
as a whole (Table 6) show that 8.8% of upland coni¬ 
fer, 4.0% of lowland conifer, 8.7% of upland hard¬ 
wood and 4.7% of lowland hardwood classes went 
into either early successional forest classes or grass¬ 
brush classes. Transitions from early successional for¬ 
est classes to mature forest classes were generally less 
than 0.5% for the ten-year period. Transitions from 
grass-brush classes to early successional conifer and 
hardwood were —5% and 15% respectively (Appen¬ 
dix 3). 

On federally managed forestland, 11.3% of the for¬ 
est area harvested between 1990 and 1995 came from 
cuts > 50 ha in size (Figure 4). This figure was 17.1% 
on private industrial forest while non-industrial pri¬ 
vate, managed state and managed county forestland 
were all similar at ~ 8.0%. Conversely, private indus¬ 
trial forest had the smallest percentage (35.6%) of re¬ 
cently cut forest area in the < 10 ha patch size classes, 
while non-industrial private forest had the greatest 
percentage (50.5%) in these smaller size classes (Fig¬ 
ure 4). 

One-way ANOVA performed on mean nearest 
neighbor (MNN) distances between early succes¬ 
sional forest patches (Figure 5) revealed a significant 
difference (P < 0.001) for one of the ownership cat¬ 
egories. Pairwise multiple comparisons showed that 
unmanaged federal forest was significantly different 
(P < 0.001) from non-industrial private forest and all 
managed ownership categories. The other significant 
difference (P = 0.011) was between managed federal 
(104.7 m ± 4.5 SE) and non-industrial private forest 
(58.6 m ± 4.2). Unmanaged federal and state forests 
had the greatest distances between patches. Within the 
unknown and managed forest categories, private in¬ 
dustrial had the greatest mean distance between early 
successional forest patches followed by federal, 
county, state and non-industrial private. 

Border lakes subsection 

An examination of the percent area of each cover type 
by patch size class indicates that strong differences 
existed between management classes in 1990 for the 
upland conifer class (Table 7). In the managed cat- 
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Table 4. Edge density (m/ha) and percent difference (%dif) 1990 to 1995 for six cover classes and 7 subsections. 



Upland conifer 

Lowland conifer 

Upland 

hardwood 


Lowland 

hardwood 

Regenerating 

forest 

Upland 

grass-brush 


Overall 


1990 

1995 %dif 

1990 

1995 

%dif 

1990 1995 

%dif 

1990 

1995 %dif 

1990 

1995 

%dif 

1990 

1995 

%dif 

1990 1995 %dif 

BL 

148.4 146.6 -1.2 

68.7 

67.8 

-1.2 

130.7 130.3 

-0.3 

9.0 

8.8 -2.7 

16.3 

24.0 

47.0 

53.6 

60.1 

12.3 

246.9 253.4 2.7 

TU 

78.8 

75.2 -4.5 

68.0 

67.4 

-1.0 

132.9 134.3 

1.0 

17.7 

17.3 -2.2 

33.7 

51.6 

53.0 

51.4 

58.2 13.4 

231.6 243.7 5.2 

LU 

117.3 

114.7 -2.2 

89.3 

89.3 

0.0 

118.5 118.2 

-0.3 

13.9 

13.5 -2.6 

15.3 

24.1 

58.0 

43.0 

50.2 16.7 

236.5 244.6 3.4 

NS 

59.0 

57.4 -2.6 

42.8 

42.3 

-1.2 

128.7 131.0 

1.8 

21.8 

21.5 -1.6 

32.3 

40.0 

23.7 

61.2 

67.6 10.5 

212.9 221.0 3.8 

NU 

83.5 

79.6 -4.7 

54.9 

54.6 

-0.5 

142.7 143.0 

0.2 

13.7 

13.4 -2.7 

41.4 

60.9 

47.2 

103.8 

113.0 

8.8 

272.5 287.3 5.4 

TL 

63.0 

60.5 -4.0 

69.8 

69.7 

-0.2 

140.3 142.1 

1.3 

19.8 

19.7 -0.8 

35.4 

50.6 

42.9 

88.8 

99.7 

12.3 

270.5 288.1 6.5 

LM 

69.4 

64.9 -6.4 

61.5 

61.4 

-0.3 

129.4 128.7 

-0.6 

24.5 

24.2 -1.5 

33.3 

52.4 

57.5 

91.3 102.0 

11.7 

257.0 273.5 6.4 


Table 5. Edge density (m/ha) and percent difference 1990 to 1995 by timber management status and year for generalized cover classes in the 
Border Lakes subsection. 



UC 

LC 

UH 

LH 

REGEN 

UG-UB 

Overall 

Managed 

1990 

127.3 

65.2 

136.2 

11.8 

26.5 

60.8 

253.3 

1995 

123.2 

63.5 

134.9 

11.4 

39.1 

71.2 

262.0 

% difference 

-3.3 

-2.6 

-1.0 

-3.5 

47.8 

17.1 

3.5 

Unmanaged 

1990 

177.1 

77.1 

124.9 

6.0 

3.4 

40.9 

249.2 

1995 

177.7 

77.0 

125.2 

6.0 

4.8 

43.0 

252.5 

% difference 

0.4 

0.0 

0.3 

-1.4 

41.7 

5.0 

1.3 

Unknown 

1990 

131.7 

56.3 

162.0 

11.2 

30.9 

80.5 

302.7 

1995 

128.0 

55.1 

160.5 

10.9 

45.8 

90.7 

313.0 

% difference 

-2.8 

-2.2 

-0.9 

-2.8 

47.9 

12.6 

3.4 


egory, 38% of the upland conifer forest was in the two 
largest size classes, (100-500 and > 500), in the un¬ 
known category 12% occurred in these size classes, 
while in the unmanaged category 72% of the upland 
conifer was in the largest size classes. Upland conifer 
area was concentrated in the smaller patch size 
classes for unknown management land, in the larger 
patch size classes on unmanaged forestland and more 
evenly distributed on managed forestland. For upland 
hardwood, managed and unmanaged showed a simi¬ 
lar distribution, with most of the area in the largest 2 
size classes, while unknown showed greatest concen¬ 
tration in the smaller classes. 


We then compared percent difference from 1990 to 
1995 for cover type area within patch size class 
within the Border Lakes Subsection and found that 
there were strong differences by management status 
in area by patch size class (Figure 6). Upland conifer 
decreased by 45 and 52% in managed and unknown 
in the > 500 ha class while unmanaged upland coni¬ 
fer showed little change across all classes. Upland 
hardwood increased in the smallest classes and de¬ 
creased in all others in the unknown management cat¬ 
egory. Upland hardwood does not occur in the > 500 
ha class in the unknown management category. In the 
managed category, upland hardwood had moderate 
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Size Class (ha) 

Figure 3. Percent change in area for patch size classes from 1990 
to 1995 for the entire study region. 

Table 6. Ten-year transition probabilities for generalized cover 
types for the entire study area. UC & LC are upland & lowland 
conifer, UH & LH are upland & lowland hardwood, GR-BR is up¬ 
land grass-brush and REGEN is forest regeneration. 


1990 

UC 

LC 

UH 

LH 

REGEN 

GR-BR 

1995 UC 

91.25 




0.13 

0.37 

LC 


95.96 



0.05 

0.21 

UH 



91.33 


0.33 

0.37 

LH 




95.34 


0.03 

REGEN 

2.36 

0.41 

5.70 

2.73 

98.09 

5.42 

GR-BR 

6.39 

3.63 

2.97 

1.93 

1.41 

93.60 


increases in most size classes, but decreased by 33% 
in the > 500 ha class (Figure 6b), which comprised 
32% of upland hardwood area in 1990 (Table 7). Up¬ 
land hardwood shows little change in the unmanaged 
category (Figure 6c). Lowland conifer had small but 
consistent decreases in all size classes (6% or less) 
where it occurs in managed and unknown categories, 


while in unmanaged, there was virtually no change 
(Figure 6a-c). 

In the early successional classes, regeneration in¬ 
creased substantially on unknown and managed land 
in all size classes except the > 100 ha classes in the 
unknown category (Figure 6d-f). Regenerating forest 
made up between 4 and 5% of the managed and un¬ 
known land in 1990 (Table 8a). On unmanaged land, 
regeneration increased by 40% in the 0-1 ha class, in 
the 1-5, 5-10, and 10-50 ha classes, increases ranged 
from 18-22% (Figure 6f), however regeneration 
makes up only 0.4% of the unmanaged land (Table 
8a). Upland grass-brush increased from 5 to 10% in 
the 3 smallest size classes in the unknown category, 
showed small decreases in the 10-50 and 50-100 ha 
classes, and increased by 25% in the 100-500 ha 
class. In the managed category, upland grass-brush 
increased in size classes where it occurred, including 
a 70% increase in the 100-500 ha class. 

Interior forest area 

Entire study area 

For the entire study area, there was a 10.5% decrease 
in interior forest area (IFA) between 1990 and 1995 
(Table 9a). On managed forestlands IFA decreased by 
12.6%. IFA on unknown management status (non-in¬ 
dustrial private) forestlands decreased by 13.5%. 
Tribal forestland had the greatest proportion of IFA 
of any management or ownership category at 53% in 
1990 and 49.5% in 1995. IFA of managed county land 
decreased in by 15.2%, second only to the industrial 
categories, which decreased by 18.8 to 24.8% (Table 
9a). 

Border lakes subsection 

Results for the Border Lakes subsection are similar 
to those of the full study area (Table 9b). Overall, the 
Border lakes subsection experienced an 8.5% de¬ 
crease in IFA from 1990 to 1995. Managed forestland 
saw a 14.3% reduction in IFA compared to a 3.8% 
decrease on unmanaged forest. Unknown manage¬ 
ment status had a 13.4% decrease in IFA. The private 
industrial forest managed category had a 25% de¬ 
crease in IFA in ~ 5 years, and it had the second 
lowest proportion of IFA in 1995 at 23.2% of mature 
forest area. In 1990, managed state forests had the 
lowest proportion of IFA of all managed forestlands 
at 26.8%. In 1995, managed private industrial forest 
had the lowest proportion of IFA (23.2%) on managed 
forestland. Tribal lands, though small in area by com- 
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H <10 ha patches 
010-50 ha patches 
□ >50 ha patches 



Federal 


State 


County Private Indust. Non-indust. Private 


Ownership 


Figure 4. Percent contribution of three patch size classes to the total forest cut between 1990 and 1995 on managed forestland and non¬ 
industrial private (unknown management status) forestlands. 
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Figure 5. Mean edge-to-edge nearest neighbor (MNN) distances between early successional forest patches by ownership and management 
status. 


parison to public ownership categories, again retained 
the greatest proportion of IFA at 64 and 59% respec¬ 
tively for 1990 and 1995. Next to tribal forestlands, 
unmanaged county forests had the highest proportion 
of IFA at 38 and 36% for 1990 and 1995 respectively. 


Managed county forestland had a 17.7% overall de¬ 
crease in IFA in ~ 5 years - second only to managed 
private industrial forest (-25.5%). Next to unman¬ 
aged private industrial forest, which represented an 
extremely small area, unknown management status 
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Table 7. Proportion of cover types represented by area class cat¬ 
egories for Border Lakes Ecological Subsection (1990). 

Border Lakes - Managed 


Area (Ha) 

UC 

LC 

UH 

LH 

REGEN 

UG-UB 

0 to 1 

11.86 

21.98 

6.24 

46.54 

29.96 

21.69 

1 to 5 

13.45 

18.44 

6.28 

25.90 

28.13 

17.85 

5 to 10 

7.18 

11.63 

4.43 

10.78 

14.66 

11.90 

10 to 50 

18.50 

29.78 

16.16 

14.27 

20.18 

32.69 

50 to 100 

11.15 

10.24 

10.06 

2.50 

4.25 

9.89 

100 to 500 

20.68 

7.94 

24.67 

0.00 

2.82 

5.97 

500 + 

17.17 

0.00 

32.14 

0.00 

0.00 

0.00 

Overall 

92,518 

42,222 

141,427 

5,563 

13,968 

40,254 

Border Lakes - Unmanaged 





0 to 1 

5.76 

32.63 

9.23 

75.87 

63.07 

48.29 

1 to 5 

5.86 

24.64 

8.98 

15.60 

21.98 

22.88 

5 to 10 

3.25 

11.94 

5.28 

3.16 

4.85 

7.73 

10 to 50 

9.36 

22.53 

17.29 

5.38 

10.11 

13.22 

50 to 100 

3.45 

2.92 

10.17 

0.00 

0.00 

2.26 

100 to 500 

12.11 

5.34 

23.98 

0.00 

0.00 

5.62 

500 + 

60.21 

0.00 

25.06 

0.00 

0.00 

0.00 

Overall 

172,785 

44,546 

124,302 

2,326 

1,391 

19,867 

Border Lakes - Unknown 





0 to 1 

18.05 

30.20 

8.16 

48.49 

31.26 

22.76 

1 to 5 

22.14 

25.12 

11.65 

34.08 

27.02 

21.90 

5 to 10 

12.61 

16.84 

10.13 

12.00 

11.43 

13.42 

10 to 50 

27.56 

24.92 

36.80 

5.43 

19.80 

26.54 

50 to 100 

7.47 

1.69 

14.26 

0.00 

4.76 

10.22 

100 to 500 

7.88 

1.24 

18.99 

0.00 

5.74 

5.16 

500 + 

4.30 

0.00 

0.00 

0.00 

0.00 

0.00 

Overall 

24,824 

9,697 

44,676 

1,607 

4,991 

15,941 


forestland (non-industrial private) had the lowest pro¬ 
portions of IFA per time-step of any other category 
(Table 9b). 


Discussion 

Landscape pattern and structure 

The increased area of early successional forest types 
across the landscape has created a more even distri¬ 
bution of cover types in this study region as indicated 
by overall increases in SWE. The increased popula¬ 
tion of early successional forest patches has also frag¬ 
mented large mature forest producing an overall de¬ 
crease in RC (Table 3). Although there are landscape 
pattern differences between the subsections and 
groups of subsections studied, they all depict a trend 


of increasing SWE and decreasing RC values from 
1990 to 1995 (Table 3). All forest management cat¬ 
egories showed overall increases in ED and SWE. 
Reductions in ED observed within mature forest 
classes were offset by large increases in ED that oc¬ 
curred within early successional types, which pro¬ 
duced an overall increase in ED for the study area 
(Tables 4 and 5). In the unmanaged forest category, 
where disturbance is largely due to wind-throw, fire 
and flooding, changes in SWE, RC and ED were 
slight compared with managed forest (Table 3). On 
managed and unknown management status lands, 
evenness of cover type distribution (SWE) is increas¬ 
ing with a corresponding decrease in forest connec¬ 
tivity (RC) and patch size. Logging operations be¬ 
tween 1990 and 1995 are primarily responsible for the 
overall effects on SWE, RC, ED and patch size mea¬ 
surements observed in this study (Plate 1). 

Upland conifer, lowland conifer, and lowland 
hardwood classes experienced decreases in ED within 
each major subsection studied, while the upland hard¬ 
wood class showed negligible to slight increases in 
ED (Table 4). This is largely due to the fact that up¬ 
land hardwoods (44% of 1990 classes) make up the 
dominant matrix of the landscape (Mladenoff et al. 
1993) when compared to upland conifer (21%), low¬ 
land conifer (15%) and lowland hardwoods (2%). 
Currently, transitions of upland hardwood to early 
successional types rarely involved the complete re¬ 
moval of individual patches from the landscape. The 
upland conifer class experienced stronger negative 
changes in ED because transitions to early succes¬ 
sional types had much greater effects on patch size 
compared to the much larger upland hardwood class. 
The relatively large decreases in ED reported for the 
upland conifers class within respective subsections is 
evident in Table 2 and in the ten-year transition ma¬ 
trix (Appendix 3) which shows that this class has ex¬ 
perienced the largest proportional transition to early 
successional types (8.75%). 

SWE and RC values indicate little change in cover 
type evenness and contagion on unmanaged forest¬ 
land in the Border Lakes Subsection. On managed 
and unknown, percent difference values were on the 
high end of those recorded for all subsections within 
the study area indicating both increased evenness and 
decreased contagion which suggests continuing frag¬ 
mentation of larger mature forest patches. The analy¬ 
sis of cover type area by patch-size class and edge 
density supports this interpretation. 
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Figure 6. Percent change in area for patch size classes from 1990 to 1995 stratified by management category for the Border Lakes subsec¬ 
tion. 


Within the Border Lakes Subsection, we found that 
patch size class distribution differed in 1990 for up¬ 
land conifer and hardwood by management status, 
indicating that these patterns were embedded on the 
landscape due to past management (Mladenoff et al. 
1993; Wallin et al. 1994). In analyzing differences in 
patch size class distribution from 1990 to 1995, we 
found that managed and unknown management status 
lands changed in a similar way to the overall land¬ 
scape, while unmanaged showed relatively little 
change across patch size class. In this subsection, 
where upland conifers are more abundant and much 
of the area is unmanaged (BWCA and Voyageurs Na¬ 
tional Park), the decrease in ED for the upland coni¬ 
fer class was substantially less than that reported for 
other subsections (Tables 4 and 7). 

We have shown that managed and unmanaged 
lands here are following different trajectories; land¬ 
scape pattern and structure in unmanaged remained 
relatively constant from 1990 to 1995, while managed 
lands showed increased fragmentation. Mladenoff et 
al. (1993) found strong differences in landscape pat¬ 


tern and structure when comparing a managed forest 
landscape to an unmanaged old-growth forest land¬ 
scape in northern Wisconsin and Michigan’s Upper 
Peninsula. Their analysis showed that on the unman¬ 
aged landscape, IFA was greater, patches were larger, 
landscape diversity was low, and landscape domi¬ 
nance high. Our comparison of landscape pattern and 
structure by management status in the Border Lakes 
showed similar results. However, the unmanaged 
lands here are not entirely undisturbed landscapes; 
portions of both the Boundary Waters Canoe Area 
Wilderness and Voyageurs National Park were man¬ 
aged for commodity production prior to receiving 
protected status (Heinselman 1996). 

Interior forest area 

Non-industrial private forestland had the lowest pro¬ 
portions of IFA of any ownership or management cat¬ 
egories (Table 9a). Three factors may be responsible 
for this trend: 1) high rates of disturbance (clear-cut- 
ting); 2) smaller average cut-units (Li et al. 1993) 
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Table 8. Area and proportion of major cover types (1990) within each timber management status for A) the Border Lakes and B) the whole 
study region. 



Managed 

Area (ha) 

% of total 

Unmanaged 

Area (ha) 

% of total 

Unknown 

Area (ha) 

% of total 

uc 

92,518 

27.54 

172,785 

47.31 

24,824 

24.40 

LC 

42,222 

12.57 

44,546 

12.20 

9,697 

9.53 

UH 

141,427 

42.10 

124,302 

34.04 

44,676 

43.91 

LH 

5,563 

1.66 

2,326 

0.64 

1,607 

1.58 

REGEN 

13,968 

4.16 

1,391 

0.38 

4,991 

4.91 

UG-UB 

40,254 

11.98 

19,867 

5.44 

15,941 

15.67 

Total 

335,953 


365,218 


101,734 



Managed 


Unmanaged 


Unknown 



Area (ha) 

% of total 

Area (ha) 

% of total 

Area (ha) 

% of total 

UC 

199,507 

17.82 

177,837 

43.80 

70,663 

11.76 

LC 

211,836 

18.93 

49,824 

12.27 

57,044 

9.50 

UH 

518,564 

46.33 

147,206 

36.25 

273,288 

45.49 

LH 

29,997 

2.68 

3,595 

0.89 

14,063 

2.34 

REGEN 

49,399 

4.41 

3,028 

0.75 

34,333 

5.71 

UG-UB 

109,996 

9.83 

24,564 

6.05 

151,385 

25.20 

Total 

1,119,298 


406,053 


600,776 



than on managed public forests (see Figures 3 and 6) 
and 3) a lack of spatially coordinated forest manage¬ 
ment practices on non-industrial private forestland in 
this region. 

Spies et al. (1994) concluded that rate of forest 
cutting has a greater effect on edge density and IFA 
than the pattern in which disturbance occurs. How¬ 
ever, other work has shown that dispersed disturbance 
patterns produce landscapes with greater edge density 
and less IFA than landscapes subjected to more ag¬ 
gregated disturbance patterns (Franklin and Forman 
1987; Li et al. 1993; Spies et al. 1994; Wallin et al. 
1994; Gustafson 1996). Timber harvest simulations 
on the Hoosier National Forest revealed that aggrega¬ 
tion of cuts in space and time produced amounts of 
edge and IFA comparable to dispersed cutting strate¬ 
gies that used half the harvest rate (Gustafson 1996) 
indicating that differences in the configuration of har¬ 
vests had greater effects on IFA than harvest rate 
(Gustafson and Crow 1996). When harvest rates are 
high and only space is considered, the advantages of 
aggregating harvest patches are soon lost (Gustafson 
and Crow 1996). 

Dense cutting patterns on non-industrial private 
forestland (Figure 5) are compounded by the fact that 
they are composed of relatively small cut-unit sizes 
(65.7% of 1990 to 1995 cut area was from patches 
^ 16 ha). In this region, square, 16 ha ( — 40 acres) 
patches are, by far, the most common parcel size con¬ 


figuration of non-industrial private ownership (Gap 
Stewardship Database, Minnesota Department of Nat¬ 
ural Resources (1998)). This is considerably smaller 
than the average managed parcel or cut-unit size (Fig¬ 
ure 6), which is consistent with what Turner et al. 
(1996) found in the Olympic Peninsula and southern 
Appalachians; that non-industrial private lands had a 
greater abundance of smaller forest patches but less 
forest cover than public lands. Using dispersed cut¬ 
ting configurations, others have found that as cut-unit 
size decreases, large interior forest patches disap¬ 
peared with increasing rapidity (Franklin and Forman 
1987; Li et al. 1993). Although a remedy for negative 
landscape effects related to spatial configuration of 
ownership may be unattainable, aggregation of small 
cut-units into larger, but fewer, cut-units would sus¬ 
tain greater IFA and less edge (Gustafson and Crow 
1994). 

Non-industrial private forest had the greatest pro¬ 
portion of its cut area in smaller patch size classes 
(Figure 4) and the most densely arranged cut-units 
(Figure 5) of unknown and managed forest categories. 
Cut-units on federally managed forestland are more 
widely spaced and have a greater proportion of the 
area in larger patch size classes than on other man¬ 
aged forest categories - except private industrial for¬ 
est where both of these patch metrics are greater (Fig¬ 
ures 4 and 5). Furthermore, IFA for each time period 
is higher on managed federal land than on unknown 
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Table 9. Interior forest area by ownership, timber management status and time for (A) the overall study region and (B) the Border Lakes 
subsection. Also, interior forest area as a proportion of total forest area per time step, and percent change in interior forest area from 1990 to 
1995. Interior forest area was calculated using a 100-meter forest edge buffer. 


A) 

Ownership 

Management 

status 

Interior forest 

1990 (ha) 

% of 1990 forest 

area 

Interior forest 

1995 (ha) 

% of 1995 forest 

area 

% diff. 1990 to 

1995 

Federal 

managed 

139,308 

32.9 

123,911 

30.7 

-11.1 


unmanaged 

109,270 

35.5 

105,163 

34.4 

-3.8 

State 

managed 

55,294 

28.3 

48,483 

26.2 

-12.3 


unmanaged 

20,866 

36.3 

19,740 

34.9 

-5.4 

County 

managed 

78,450 

27.1 

66,510 

24.2 

-15.2 


unmanaged 

1,235 

39.2 

1,164 

37.4 

-5.7 

Private Industrial 

managed 

10,802 

28.6 

8,772 

25.3 

-18.8 


unmanaged 

2,118 

21.3 

1,593 

17.7 

-24.8 

Tribal 

managed 

8,746 

53.0 

8,181 

49.5 

-6.5 

Non-industrial 

Private 

unknown 

80,579 

19.4 

69,729 

17.6 

-13.5 

All Owners 

managed 

291,807 

30.3 

255,126 

27.9 

-12.6 


unmanaged 

133,488 

35.3 

127,660 

34.1 

-4.4 


unknown 

80,579 

19.4 

69,729 

17.6 

-13.5 

B) 

All combined 

505,874 

28.8 

452,515 

26.9 

-10.5 

Federal 

managed 

46,288 

28.9 

40,562 

26.9 

-12.4 


unmanaged 

107,947 

35.6 

103,884 

34.5 

-3.8 

State 

managed 

18,872 

26.8 

15,933 

24.2 

-15.6 


unmanaged 

13,834 

37.6 

13,285 

36.5 

-4.0 

County 

managed 

11,983 

29.5 

9,860 

25.9 

-17.7 


unmanaged 

1,180 

38.0 

1,112 

36.1 

-5.8 

Private Industrial 

managed 

2,302 

27.4 

1,715 

23.2 

-25.5 


unmanaged 

47.3 

8.7 

46.8 

8.6 

-1.1 

Tribal 

managed 

2,155 

64.0 

1,907 

59.1 

-11.5 

Non-industrial 

Private 

unknown 

19,901 

24.6 

17,237 

22.5 

-13.4 

All Owners 

managed 

81,272 

28.8 

69,669 

26.3 

-14.3 


unmanaged 

123,009 

35.8 

118,327 

34.7 

-3.8 


unknown 

19,901 

24.6 

17,237 

22.5 

-13.4 


All combined 

224,181 

31.7 

205,232 

30.0 

-8.5 


and other managed forestland, which likely reflects 
ownership-related differences in management (Turner 
et al. 1996). With the exception of managed private 
industrial forest, the rates of disturbance on managed 
forestland are similar across all ownership categories 
(Table 1). Interestingly, the disturbance rate between 
managed federal and non-industrial private forest is 
virtually identical at 0.95% per year, while IFA of the 
latter was the lowest of all ownership categories for 
both time periods. Lastly, the variability of cut-unit 
spacing alone may have an effect on IFA as Gustafson 
(1996) simulations suggested. Private industrial for¬ 
est had greater variability in cut-unit spacing than 


managed federal forest (Figure 5), while cut-unit size 
remained similar. The proportion of IFA on private 
industrial forestland (25.3%) was similar to that on 
managed federal land (30.7%) even though harvest 
rate was —43% greater on the former suggesting a 
cut-unit pattern effect on IFA. 

IFA is strongly related to disturbance patch size. 
Large matrix forest patches (> 500 ha) have an im¬ 
portant influence on landscape pattern and structure; 
they maximize forest interior (habitat) and connectiv¬ 
ity (Mladenoff et al. (1993); Crow et al. 1999). Our 
analysis of patch size distribution, IFA, RC and MNN 
shows an overall trend toward less IFA and decreased 
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connectivity across the managed forest landscape in 
Northeastern Minnesota. 

Influence of management status and physical 
environment on landscape pattern 
Soil and landform characteristics have been shown to 
influence forest landscape pattern and structure in 
Northern Lake States forests (Pastor and Broschart 
(1990); Mladenoff et al. (1993); Crow et al. 1999). 
Crow et al. (1999) also demonstrated that the interac¬ 
tion of management differences related to ownership 
and soil-landform ecoregions influence landscape pat¬ 
tern. Our analysis suggests that although there is vari¬ 
ation in landscape metrics (SWE, RC, ED) related to 
ecological subsection areas, management has a pro¬ 
found influence on landscape pattern and structure in 
the study area. The same shift toward increasing 
cover type evenness (SWE), decreasing cover type 
contagion (RC) and increasing edge density (ED) (es¬ 
pecially in early successional types) occurred across 
all subsections. Within the Border Lakes Subsection, 
analyses (SWE, RC, ED, patch size) were stratified 
by management status and indicate pattern and struc¬ 
ture differed between managed and unmanaged at the 
onset of the study. Results from 1990 to 1995 indi¬ 
cate continued divergence with managed and un¬ 
known showing the same shifts toward increasing 
fragmentation that occurred over the whole study 
area. Although we did not measure IFA by subsection, 
IFA decreased substantially in all managed forest cat¬ 
egories, which span all subsections in the study area. 
IFA also varied by ownership within the managed 
forest category. These results suggest that the patterns 
imposed by management are to some degree overrid¬ 
ing landform control and creating a more homoge¬ 
neous landscape pattern across the study area. Turner 
et al. (1996) determined that ownership had a strong 
influence on both landscape pattern and change in the 
Olympic Peninsula and Southern Appalachian High¬ 
lands study regions. They also noted that land own¬ 
ership may have important effects on ecological pro¬ 
cesses because of the strong influence of ownership 
on landscape pattern. Our results in northeast Minne¬ 
sota also show that management, as affected by own¬ 
ership strongly influences landscape pattern. Turner et 
al. (1996) suggest that ownership should be explicitly 
considered when the future health of ecosystems is 
assessed. 


Biological diversity and forest ecosystem health 

Our results show that forest patch size, interior forest 
area and connectivity have decreased with a concom¬ 
itant increase in edge density. Additionally, current 
management practices may favor hardwood regenera¬ 
tion over conifer species. These forest landscape 
trends may have important implications for biologi¬ 
cal diversity and forest ecosystem health. 

Landscapes with small patches, high amounts of 
edge and significant areas in early successional cover 
types likely provide benefits to valued wildlife spe¬ 
cies such as ruffed grouse and white tailed deer 
(Hunter 1990). However, management for high edge 
density and early successional cover types may have 
significant effects on regional biological diversity. 
The current high population of white-tailed deer and 
their browsing on eastern white pine and northern 
white cedar is the single most important factor in the 
current low regeneration rates of these two tree spe¬ 
cies (Anderson et al. 2002). The research results on 
the effects of edge-interior relationships on bird spe¬ 
cies populations are inconclusive for primarily for¬ 
ested landscapes. However, at the regional scale, in¬ 
creasing proportions of forest edge may compromise 
the viability of some edge-sensitive forest birds (Tem¬ 
ple and Flaspohler 1998). 

Changes in landscape pattern and structure may 
also affect forest health in this region. Recent work 
has shown that jack pine budworm (Kouki et al. 
1997) and forest tent caterpillar (Rolands and Kaup 
1995) outbreaks are more severe in highly patchy for¬ 
ested landscapes with high edge density and small 
patch sizes. 

Our results suggest that current and past forest 
management systems may have adverse effects on bi¬ 
ological diversity and overall forest ecosystem health 
in this region. However, the legacy of past forest 
management practices may be difficult to erase with¬ 
out significant reductions in disturbance rates and/or 
changes in minimum harvest age (Wallin et al. 1994). 

Usefulness of remote sensing for monitoring 
landscape pattern and structure 

Landsat TM data has proven to be a useful tool for 
mapping regional forest cover and for detecting 
changes in forest cover through time. Image differ¬ 
encing techniques using the normalized difference 
SWIR/NIR ratio was an efficient way to detect spec- 
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tral changes resulting from alterations in forest cano¬ 
pies. 

In the past, landscape pattern and structure mea¬ 
surements for large landscapes were performed using 
coarse land cover data from one time period (Mlad- 
enoff et al. 1997). This one-time snapshot using 
coarse resolution data is of limited use for determin¬ 
ing regional landscape management needs. Others 
have focused on detailed, photo-based land use 
change studies covering multiple time-steps for rela¬ 
tively small regions (White and Mladenoff 1994). Al¬ 
though these studies produce more precise pictures of 
landscape scale trends through time, it is uneconomi¬ 
cal to perform such work on the scale necessary for 
many resource management decisions. 

Satellite data affords both appropriate scale and 
time line for regional studies of landscape pattern and 
structure. One advantage to the use of remote sensing 
data for these inventory purposes is that one can ac¬ 
curately determine forest type as well as area affected 
and do so on a regular basis. Future effort is needed 
to develop a more detailed remote sensing-based time 
line of forest management in northern Minnesota and 
elsewhere. Annual observations using TM data could 
facilitate analysis of forest harvesting behavior in re¬ 
lation to variable market values of wood (Turner et 
al. 1996). Within the time frame we studied, the mar¬ 
ket value of aspen in this region rose from USD 7.00/ 
cord to roughly USD 21.00/cord, and pulpwood value 
for many other tree species were at all time highs 
(Minnesota Forest Resources Council 1999). By ac¬ 
quiring TM data from before and after the time inter¬ 
val studied here, we can build a more comprehensive, 
species-by-species account of what has taken place in 
Minnesota’s forests. Remote sensing technology will 
allow us to make more informed decisions regarding 
landscape-based forest resource planning and policies 
that affect both forests and the wildlife dependant 
upon them. 

Lastly, it should be noted that conifer regeneration 
appears to be under represented in the ten-year tran¬ 
sition matrix (Appendix 3). Based on field observa¬ 
tions and 1984 TM data, we believe most of the 
younger conifer regeneration patches that originated 
in the late 1970s and early 1980s, as a result of cut¬ 
ting 1933-34 origin conifer plantations, are seen as 
either grass or brush using 1990 TM data. Conifer re¬ 
generation (< 15 years old) does not form closed can¬ 
opy patches like hardwood regeneration patches do, 
which allows brush and grass in the young conifer 
patch matrix to dominant the TM spectral signature. 


When conifer patches approach crown closure, their 
TM spectral signature becomes more distinct, which 
may explain why a substantial amount of the 1990 
grass and brush class moved to conifer regeneration 
by 1995. 

Summary 

We used remote sensing technology to classify forest 
cover to a near-species level for a 2.76 x 10 6 ha area 
and, by using image ratio differencing techniques us¬ 
ing the normalized difference SWIR/NIR ratio 
(NDSN), we were able to efficiently identify changes 
that occurred in the forest cover during a 5-year pe¬ 
riod. 

Overall, 4.2% of the mature forested area was 
changed in the five year period studied. Of this, 5.3% 
of the spruce-fir, 4.7% of the aspen-birch, 4.6% of the 
jack pine and 3.2% of the black spruce cover type 
changed to early successional types. Forest harvest¬ 
ing on managed forests was highest on private indus¬ 
trial land at —1.7% per year and lowest on tribal 
land, —0.55% per year. Managed State forests had 
the greatest rates of disturbance of all publicly man¬ 
aged forestland (— 1.1% per year). 

We stratified the landscape by ecological subsec¬ 
tion, ownership and timber management status and 
found that landscape structure varied significantly by 
each of the three stratification categories. The most 
dramatic landscape pattern effects appear to be driven 
mostly by management status and ownership. Tribal 
forestlands had the highest percent of IFA, while non¬ 
industrial private forests had the lowest. Non-indus¬ 
trial private lands have a greater proportion of their 
cut area in smaller patch sizes than managed public 
or private industrial forests. Differences in IFA ob¬ 
served between managed public, non-industrial pri¬ 
vate and private industrial forests are most likely due 
to differences in cut-unit size, harvest spacing pat¬ 
terns and, to a lesser degree, rate of harvest, and 
hence, reflects the management and land-use strategy 
of the owner (Turner et al. 1996). 

Cover type area by patch size class analyses 
showed a general trend of decreasing area in most 
patch size classes for upland conifer, lowland conifer 
and lowland hardwood. Upland conifer showed its 
greatest decrease (8%) in the largest size class (> 500 
ha). Upland hardwood showed an increase in all patch 
size classes except the largest class, which showed a 
13% decrease in area. 
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In Summary: 

• Harvest rates vary by timber management status, 
but show less variability by ownership category. 

• Fragmentation is greatest on non-industrial pri¬ 
vate forests. Overall fragmentation trends may 
influence both biological diversity and ecosys¬ 
tem functions. 

• Patterns of ownership affect landscape patterns 
(non-industrial private versus managed public 
forests). 

• Our analysis of RC, patch size distribution, IFA, 
and MNN shows an overall trend toward less 
IFA and decreased connectivity across the man¬ 
aged forest landscape in Northeastern Minne¬ 
sota. 

• Cut-unit size and arrangement appear to have a 
greater effect on IFA than observed rate of har¬ 
vest. 
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Appendix 1. Confusion matrix for the 1990 forest classification (WRS 27/27) 


0.739 Overall Accuracy = 75 % 
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Reference Data 
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4 
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0.92 

0.95 
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Overall accuracy of change classes 

88.8 


KHAT 

0.871 






Appendix 2. Accuracy assessment for 1990 to 1995 change classes 
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Appendix 3. Ten-year transition probabilities for the whole study region. Diagonal elements are retention frequencies and off-diagonal ele¬ 
ments are transition frequencies. 
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